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Abstract 

We present a framework to analyze effects of new physics beyond the standard 

model on B-B mixing and CP violation in B decays in a model- independent 

manner. Assuming that tree level decay amplitudes are dominated by the 

standard model ones, new physics contribution to the B-B mixing can be 

extracted from several measurements at B factories. Using this framework, we 

show the present constraint on new physics contribution to the B-B mixing, 

and illustrate constraints expected to be given by future experiments at B 

factories. We also point out a possibility that CP asymmetries in B — ► ipKs, 

B — ► 7T7T, and B — ► DK modes look consistent with the standard model, even 

if a large new physics contribution is present in the B-B mixing. 
12.90. +b, 13.25.Hw, 14.40.Nd 
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Physics of B meson provides several tests of the standard model and could give insights 
into new physics beyond it. Especially, in the standard model, the test of unitarity of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix |I| is most important. As shown in Fig. [1], the 
unitarity of the CKM matrix is graphically expressed by a triangle. In the standard model, 
the lengths of the sides are related to several decay rates and/or the magnitude of the B-B 
mixing, while the angles are related to several CP asymmetries. When these quantities are 
measured at future B factories, we will be able to test the unitarity of the CKM matrix by 
looking whether the triangle is closed or not. 

Quantitatively, the analysis will be done along the y 2 test: All the observables such as 

decay rates and CP asymmetries are represented by the standard model parameters like the 

CKM parameters and some hadronic parameters like the B meson decay constant The 

X 2 is calculated with experimental results of these quantities and knowledge of the hadronic 

parameters as a function of the standard model parameters. Then, this x 2 is minimized 

by varying the standard model parameters, and we can see whether the standard model is 

consistent or not depending on the obtained y?L- . 

mm 

Although it is quite straightforward and powerful, this method has some defects. Even 
if a larger value of yL™ * s f° un d> this method itself does not tell us anything other than 
that the standard model is doubtful. No quantitative information about new physics can 
be obtained. Moreover, even if the standard model seems to be consistent, there may be 
new physics which evades the above y 2 test. From these points of view, it is desirable to 
introduce some new physics effects into the analysis and to see how they are restricted by 
experiments. 

For this purpose, two approaches are possible: One is a mo del- dependent approach, in 
which we specify a model of new physics and analyze all data in it. Although each model 
needs a specific analysis, this approach has a predictive power if the model does not have 
too many parameters. Another approach is a model-independent approach, in which we 
do not introduce any model. The data are analyzed based on rather general principles 
or assumptions. If we can parametrize effects of new physics by a few model-independent 
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quantities, this approach is quite suitable to select appropriate new physics among many 
possibilities based on experimental information. We pursue the latter approach in this paper. 

Weak interactions of B meson are described by |AB| = 1 and |A5| = 2 amplitudes, i.e. 
B decay and B-B mixing amplitudes respectively. In the standard model, the |A£?| = 1 
processes occur through the tree and penguin diagrams at the quark level. Effects of new 
physics tend to appear in a penguin diagram because it is a loop diagram. In fact, in some 
classes of new physics like SUSY, tree-level B decay amplitudes are hardly affected by new 
physics. On the other hand, it is difficult to exclude effects of new physics in the penguin 
diagrams. 

In order for the model-independent approach, we classify B decay processes into two 
classes. Although we are considering general cases, it is useful to introduce class I and 
class II based on the properties of the quark-level amplitudes in the standard model. The 
class I processes do not have the penguin part in its quark-level amplitude, and the class II 
processes do. The processes of b — > qdq u q' u tyP e ) where q<i represents a down-type quark, and 
q u and q' u denote different up-type quarks, are class I processes. The ordinary semileptonic 
decays are also regarded as class I processes because their amplitudes have no penguin part. 
While, the processes of b — > q d qq type, where q denotes a generic type of quark, are class II 
processes. If q is an up-type quark, b — > qdqq consists of the penguin part and the tree part. 
If q is a down-type quark, b — > qdqq contains only the penguin partQ. 

In the following, we assume that the class I processes are described by the standard 
W-exchange diagrams. In addition, we assume that the unitarity of the CKM matrix is 
saturated by the first three generations. This assumption is necessary to determine the 
quark-mixing parameters from class I processes and to apply them to evaluate the standard 

x The processes b — ► dds and b — ► ssd do not appear in the standard model within the lowest order 
of the weak interaction. Even if these processes are caused by new physics, they do not affect the 
following discussions. 
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model contribution to the B-B mixing as is described below. 

New physics can also contribute to the B-B mixing amplitude M 12 . We can write it as 

M 12 = \M 12 \ e^ M = Mf 2 M + Mf 2 EW , ( 1 ) 

where M^ 2 M is the standard model contribution and M^ EW represents contributions from new 
physics. Mf 2 EW is the model-independent parameter which describes effects of new physics 



in the B-B mixing;. Note that Eq. (I) 

is enough to determine a phase convention for M 



NEW 
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once the phase convention in evaluating Mf 2 M is fixed, because the relative phase between 
Mf 2 M and M^ EW is a physically meaningful quantity. We use the following expression of 



Mf 2 M : 



M X S 2 M = |^ m > B /^^(W) 2 5KR), (2) 



where 



(V; d V tb ) 2 = X G A 2 (l-p + iV ) 2 , (3) 
in the Wolfenstein parametrization M, and 



S(x) = x 



1 9 3 

- + 



x 3 



2(1 -x) 



lnx, (4) 



4 4(1 -x) 2{l-x) 2 _ 
is the Inami-Lim function ||. 

In the following, we see how M^ EW is constrained by the present experiments, and 
how the constraint will be improved by future experiments. For this purpose, we have to 
determine the Wolfenstein parameters A, p, and rj simultaneously with Mj^ EW . The ordinary 
analysis in which the standard model is assumed cannot be applied if we consider effects of 
new physics. However, according to the above two assumptions, i.e. the dominance of the 
standard model contributions in the class I processes and the three-generation unitarity, we 
can extract information about the Wolfenstein parameters without being bothered by the 
unknown new physics in the class II processes which is not parametrized in our analysis. 

First, let us consider the semileptonic decay of the B meson, B —* X c £u, which is one of 
the class I processes and free from new physics. We can determine the CKM matrix element 



\Vcb\ = X 2 A from its width. Secondly, the charmless semileptonic decay of the B meson, 
B — > X u iv is also considered to be free from new physics, and we can obtain a constraint 
on |V^b/\46| = A a/ p 2 + rj 2 from its rate. These exhaust the presently available constraints 
which are not affected by new physics in our framework. 

The remaining constraint relevant to our analysis is the one given by the observation of 
B-B mixing. From the experiments, we obtain a constraint on the absolute value of the 
B-B mixing amplitude, IM12I = |Mf 2 M + M£j, EW |. Since the possible range of Mf 2 M is limited 
by the above mentioned constraints on the Wolfenstein parameters and the knowledge of 
the top quark mass and the hadronic matrix element, we can obtain an allowed region of 
^new ^ no ^. j nc i uc i e information from the K-K mixing into our analysis because 
we need to introduce another model-independent parameter, M^ EW for the K-K mixing. 

Here, we present a result of the analysis which follows the above strategy. We summarize 
the inputs in Table |. A is determined by the semileptonic kaon and hyperon decays [|J, 
which are free from new physics in our framework. A and V p 2 + rj 1 are constrained by the 
semileptonic B decays as mentioned above ||. IM12I is obtained from the measurement of the 
mass difference Am || . The top quark mass mt is given by the direct observations [|7| . 
IbVBbVb is estimated by several theoretical methods. We take the value used in Ref. [|J. 
We omitted experimental errors and/or theoretical uncertainties of A, A, mt, and Jb^/Bbtib 
for simplicity. 

Fig. |2| shows the experimentally allowed region and the standard model prediction for 
Myi in the complex M 12 plane. The circles whose centers are at the origin show the direct 
experimental constraint on \Mxi\ in Table |, and the distorted circles show the standard 
model "prediction" which is calculated by Eq. (0) with the inputs in Table |I[ A vector 
from a point in the region of the standard model prediction to a point in the experimentally 
allowed region corresponds to a possible complex value of . From this figure we 

can see that the value of M^ EW is certainly limited. By integrating out p and rj in the 
maximum likelihood method, we can obtain an allowed region of M^ EW . The result is 
shown in Fig. |||(a), where we show the line that x 2 = Xmm + 4.61, which is usually regarded 
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as 90% confidence level. This figure tells us that the contribution of M^ EW to the total B-B 
mixing amplitude M 12 can be comparable with or larger than that of iWf 2 M , besides that the 
standard model is completely consistent. Also note that a positive value of ReM^ EW tends 
to be disfavored. 

Now, let us discuss effects of new physics on CP asymmetries in B decays which will be 
measured in future experiments at B factories. We start with the class I process B — > DX S , 
for which we can use the standard model decay amplitudes. For the charged B meson mode 
B — > DXf, since this decay mode is free from new physics provided that no significant 
CP violation in D°-D° mixing is caused by new physics, the method by Gronau and Wyler 
H can be applied to extract information of the angle 3 defined in Fig. |l|. We can uniquely 
determine sin 2 3 using two strange states which make different phase shifts from each other, 
e.g. K and K*, for X s . For the neutral B meson mode B°(B°) — > DK$ [HJ, the situation 



is a little different owing to the new physics contribution to the B-B mixing, M^ EW . From 
the time-dependent decay rate, however, we can uniquely determine sin0 M , sin(0 M + 20 3 ), 
and sin 2 3 Q. 

For B — »■ 7T7T process, the decay amplitude can be affected by new physics because it 
is a class II process. However, we can extract the AI = 3/2 part of the amplitude by the 
isospin analysis fll |. Since the AI = 3/2 part of the amplitude does not contain the penguin 



contribution, it can be regarded as a class I process. Although this analysis depends on the 
assumption that the penguin part of the amplitude is AI = 1/2, it seems to be plausible even 
in the presence of new physics. In other words, we do not expect that the electromagnetic 
penguin, the electroweak penguin, and the box diagram play important roles in this process 
even in the presence of new physics as in the standard model [0. In such a case, we can 
extract sin(0Af + 20 3 ) with fourfold ambiguity. If the penguin contribution turns out to be 



2 In the standard model, which is considered in Ref. [ 10 1 , 4>m defined in Eq. (||) corresponds to 
20i in our phase convention. 
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small, this fourfold ambiguity disappears. We can also resolve it by comparing with the 
sin(0 M + 20 3 ) obtained from B° (B°) -> DK S . 

The last process discussed here is B°(B°) — ► ipKs, which is classified as a class II 
process. In the standard model, the penguin contribution in this process does not cause 
any problem, because the weak phases of the tree and the penguin contributions are the 
same. However, they could differ from each other in the presence of new physics, so that 
we cannot extract information about 0m from this process in general. Nevertheless, it is 
desirable to include this mode in our analysis because it is expected to be precisely measured. 
In the following, we assume that there is no significant penguin contribution which has a 
different CP- violating phase from that of the standard model in this process. With this 
assumption, we can uniquely determine sin0 M from the time-dependent decay rate. Note 
that the above assumption can be checked by comparing the obtained sin0M in this process 
with that obtained in B° (B°) -> DK S process. Also, it can be tested by looking for the 
CP-violating rate differences (the direct CP violation) which may be seen in B ± — > ipK ± 
and cos(Ami) term in the time- dependent rate of B°(B°) — > ipK$ decay depending on the 
relevant phase shifts. Table [TT] summarizes the above arguments on the several modes of 
studying CP violation. 

Now, let us illustrate how measurements of CP asymmetries in these modes at B facto- 
ries constrain the model- independent parameter M^ EW . The inputs for this illustration in 
addition to those in Table | are given in Table |T|. The central values of these inputs are 
calculated by putting M^ EW = and (p, rf) = (0.23,0.29) which is a typical point allowed 
in the standard model analysis ||. The errors in this table are taken from Ref. [13| . 



At an earlier stage of a B factory run, we will observe only the CP violation in B°(B°) — > 
ijjKg mode. In Fig. 0(b), we show the expected constraint on M^ EW from this mode which 
gives the constraint on sin 0^ and the information given in Table |I[ The expected constraint 
is not much stronger than that in Fig. 0(a), because 03 still remains free. In Fig. 0(c), we show 
the expected constraint by adding the information from B^ — > DXf i.e. the constraints on 
sin 2 03. In this case, we have eight solutions for M^ EW owing to twofold ambiguity in 0ju 



determination from sin 0m and fourfold ambiguity in 03 determination from sin 2 03. These 
solutions are indicated by the dots in Fig. 0(c). The rather larger allowed regions in Fig. 0(c) 
consist of those surrounding these eight solutions. Fig. 0(d) shows the expected constraint 
from all the inputs given in Table | and fTJ. We have four solutions even in this case. If 
allowed regions obtained with real data do not contain the origin (M^ EW = 0), we can 
not only conclude that the standard model is excluded but also determine the new physics 
contribution quantitatively. 

An important observation from this illustration is that we cannot exclude the possibil- 
ity of significantly large M^ EW even if all the measurements considered here seem to be 
consistent with the standard model. 

We have neglected the uncertainties of A, m t , and /bVBbVb in the above analysis. 
The allowed regions in Figs. 0(a)— (d) are changed if we vary these values. However, we 
expect that these uncertainties will be reduced enough by future experimental and theoretical 
developments. 

In conclusion, we presented a framework of analysing the B-B mixing and the CP vio- 
lations in B decays in a model-independent manner. We introduced the model-independent 
parameter M^ EW and showed the constraint on it from the presently available experimental 
data. We also illustrated how M^ EW would be constrained by the future experiments of CP 
violation at B factories. We found that there remain some non-trivial solutions of M^ EW 
even for the inputs corresponding to the standard model case (Mj\j, EW = 0). This means 
that we cannot exclude the possibility of significant contribution to the B-B mixing from 
new physics even if CP asymmetries in all the modes of B — > ipKs, B — > DX S , and B — > tttt 
seem to be consistent with the standard model. 

The authors would like to thank C. S. Lim for useful discussions. 
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TABLES 



A 


0.220 @ 


\V cb \ = X 2 A 


0.038 g 


m t 


174 GeV § 




165 MeV §§ 


\V ub /V cb \ = \^fP + V 2 


0.08 ± 0.02 § 


Am = 2\M 12 \ 


0.462 ± 0.026 ps" 1 § 


TABLE I. The inputs 


corresponding to the present experimental and theoretical knowledge. 


B ± -> DXf 


sin 2 03 


B°(B°) -» DK S 


sin0 M , sin(0 M + 20 3 ), sin 2 3 


B — > TTTT 


sin(0Af + 2^3) (fourfold in general) 


B°(B°) -» ^ 


sin 0m (assuming no penguin which has a non-standard phase) 


TABLE II. Quantities obtained from several processes of CP violation search. 


sin (j) M 


0.66 ±0.08 


sin 2 03 


0.62 ±0.25 


sin(0M + 203 ) 


0.57 ±0.17 



TABLE III. The inputs used for the illustration of the expected constraints from CP violation 
experiments at B factories. 
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FIGURES 



FIG. 1. The unitarity triangle. 

FIG. 2. The experimentally allowed region and the standard model prediction for M\2- The 
circles whose centers are at the origin show the direct experimental constraint, and the distorted 
circles shows the standard model prediction. The dashed lines show the standard model predictions 
for several fixed 4>s values. (The angle <f>3 is defined in Fig. 

FIG. 3. The allowed regions of M^ EW . The lines that \ 2 = Xmin + 4 - 61 are shown, (a) The 
presently allowed region corresponding to the inputs in Table |. (b) The expected constraint from 
B — > ipK$ in addition to the inputs in Table |. (c) The expected constraint from B — ► ipKs and 
B^ — > DXf 1 in addition to the inputs in Table Q. The dots represent the solutions obtained from 
the central values of the inputs, (d) The expected constraint from all the inputs listed in Table | 



and III. The dots mean the same as (c). 



11 




Fig. J 
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